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Abstract 

The evolution of a Universe confined onto a 3-brane embedded in a 
I ' five-dimensional space-time is investigated where the cosmological fluid 

on the brane is modeled by the van der Waals equation of state. It is 

■ shown that the Universe on the brane evolves in such a manner that three 
\ distinct periods concerning its acceleration field are attained: (a) an initial 

f*"^ . accelerated epoch where the van der Waals fluid behaves like a scalar field 

t-H ■ with a negative pressure; (b) a past decelerated period which has two 

, contributions, one of them is related to the van der Waals fluid which 

behaves like a matter field with a positive pressure, whereas the other 
contribution comes from a term of the Friedmann equation on the brane 
^) , which is inversely proportional to the scale factor to the fourth power and 

can be interpreted as a radiation field, and (c) a present accelerated phase 
»' | | due to a cosmological constant on the brane. 

Key words: brane cosmology, van der Waals fluid, acceleration field 

■ According to the cosmological observations one can distinguish three dis- 
\ tinct periods for the Universe that are related to its acceleration field. The 

5^ ■ first period refers to an accelerated epoch dominated by a scalar field where a 

rapid expansion of the Universe characterizes its inflationary phase. The next 
one is related to a decelerated phase dominated by matter fields. This period 
is followed by a return to an accelerated epoch dominated by a cosmological 
constant or dark energy. 

Recently several authors have investigated the evolution of the Universe 
within the framework of the one-brane model of Randall and Sundrum ;1. where 
the Universe is confined onto a 3-brane, which is a hyper-surface embedded in a 
five-dimensional space-time called bulk. For recent reviews on this subject one 
is referred to the works of Brax and van de Bruck 2 and Langlois [3| and the 
references therein. 

In the present work we investigate the evolution of a Universe confined onto 
the 3-brane embedded in a five-dimensional space-time in order to describe 
the three distinct periods of the Universe beginning with an accelerated phase 
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passing through a decelerated epoch and returning to an accelerated period. 
For that end we model the cosmological fluid on the brane by the van der 
Waals equation of state. The use of the van der Waals equation of state in 
cosmological problems was first proposed by Capozziello and co-workers 
who recognize that this equation of state could describe the transition from a 
scalar field dominated period to a matter dominated epoch without the need 
of introducing scalar fields. Recently a model for the Universe as a mixture of 
a van der Waals fluid with dark energy (modeled as quintessence or Chaplygin 
gas) was proposed in the work [S] - within the framework of a four-dimensional 
space-time theory - in order to describe the transition from the accelerated- 
decelerated-accelerated periods of the Universe. 

In this work we show that the Universe on the brane modeled by the van 
der Waals equation of state evolves in such a manner that the three distinct 
periods concerning its acceleration field are attained. The initial accelerated 
epoch is due to the van der Waals fluid which behaves like a scalar field with 
a negative pressure. The past decelerated period has two contributions: one 
of them is related to the van der Waals fluid which behaves like a matter field 
with a positive pressure, whereas the other comes from a term of the Friedmann 
equation on the brane which is inversely proportional to the scale factor to the 
fourth power and can be interpreted as a radiation field. The present accelerated 
period is due to a cosmological constant on the brane. 

In the one-brane model of Randall and Sundrum the Universe is confined 
onto a hyper-surface - called brane - which is embedded in a five-dimensional 
space-time with coordinates (t, x 1 , x 2 , x 3 , y), called bulk. The brane is located 
at y = and the line element which describes a spatially flat, homogeneous and 
isotropic Universe on the brane is given by (see, for example, the reviews [21 El) 

ds 2 = gMN<ix M dx N = n(t, y) 2 dt 2 — a(t, y) 2 5ijdx i dx : ' — dy 2 , (1) 

where gMN is the metric tensor with signature (+,—,—,—,—). The two func- 
tions n(t,y) and a(t,y) are determined from the five-dimensional Einstein field 
equations that read 

Rmn — -RguN + — —k 2 T M n- (2) 

Above, Rmn is the five-dimensional Ricci tensor, R = R M m its trace, A denotes 
a bulk cosmological constant, k 2 = 8irG5 is related to the five-dimensional 
gravitational constant G5, and Tmn is the five-dimensional energy-momentum 
tensor. 

The five-dimensional energy-momentum tensor Tmn is decomposed into a 
sum of two terms: one refers to the bulk whereas the other is related to the 
brane, i.e., 

rpM rpM I I rpM I CQ\ 

J N = 1 AM bulk + J JVIbrane- \o) 

The energy-momentum tensor in the bulk is given by 

r M jv|bulk = diag( ( o B ,-p_B,-p_B,-p B ,p B ), (4) 
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where the energy density ps and the pressure in the bulk do not depend 
on the y coordinate. Furthermore, the energy-momentum tensor on the brane 
reads 

TMjvlbranc = S(y)(ag flL , + T^S^mS^n- (5) 

On the brane there exist two contributions: one is related to the constant tension 
on the brane ct, whereas the other refers to the energy-momentum tensor of the 
cosmological fluid which is written as 



T^v = diag 0, -p, -p, -p). 



(G) 



Above, p and p denote the energy density and the pressure of the cosmological 
fluid on the brane, respectively. 

For the line element given by Q , the components of the Einstein field equa- 
tions © become 
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A = k 2 T v . 



(8) 
(9) 

(10) 



In the above equations the dot and the prime refer to a differentiation with 
respect to the coordinates t and y, respectively. 

We follow the work [7] and introduce the abbreviations at y = 0: ao(t) = 
a(t,y = 0) and no(t) = n(t,y = 0). Hence, on the brane the line element is 
written as 

ds 2 = n (t) 2 dt 2 - a {t) 2 8 lj dx t dx J . (11) 

By imposing the gauge condition n (£) = 1, the time coordinate becomes the 
proper time on the brane and ao(t) is identified with the scale factor. 

If one assumes that there exists no matter flow along the fifth dimension 
(see, e.g. |S])- so that Tq v = - the integration of equation leads to 



n(y,t) 



d (t) 



(12) 



We integrate equations 10 and (JHJ) over — e < y < +e, and by taking the 
limit e — > 0, it follows the so-called junction conditions: 



a 



a + —p + n z p, 



(13) 
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where [[/]] = /(0+) — /(O— ) denotes the jump of the function / across the brane. 

If we take into account in equation 10 the junction conditions (|13(l together 
with the condition of no matter flow along the fifth dimension - i.e., Tq v = 
we get the conservation equation for the energy density on the brane 

p + 3^(p+p) = Q. (14) 
a 

In order to analyze the system of Einstein field equations 10 - (|10|l in the 
bulk, we follow the work [H] and introduce the function 

F(t 7 y)^(a'a) 2 ~^. (15) 
In terms of the function F(t, y) the 00, yy and ij - components in the bulk read: 

F' = -^(A + ^), F = -^f(A + K 2 p B ), (16) 



F' 



2aa\K- K 2 p B ), (17) 



respectively. Now we differentiate (tTB|) i with respect to time t and (|TT)|l 9 with 
respect to the coordinate y, and get that ps does not depend on time t. Fur- 
thermore, from the differentiation with respect to time of the expression F' /a' 
obtained from Ijl6|) i it follows an equation when compared with Ijl7|l implies 
that pb = —Pb- Finally, the integration of Ijl6|h with respect to the coordinate 
y, leads to 

F = ( a >a) 2 -(™f = -^(A + K 2 p B ) + C. (18) 
n z 6 

Above, C is a constant since ps is time independent. 

The Friedmann equation on the brane is obtained from i|18|) by considering 
the limit y — > and the junction conditions H13|) . yielding 

a J 18 36 36 

From the Friedmann equation one can determine the acceleration equation 
on the brane. Indeed, the differentiation of (|19fl with respect to time, yields 

- - ~*(P + 3p) - ^P(2p + 3p) + ^y° 2 + 6A + 6k 2 Pb ) - - A , (20) 
ao 36 36 36 

thanks to the equation (|14fl for the energy density on the brane. This equation 
can be obtained also from equation l|10l) by using the junction conditions (|13f) . 

Now we have a system of differential equations for the determination of 
the energy density p(t) and of the scale factor ao(t) on the brane which is 
composed by the conservation equation for the energy density 1141) and by the 
Friedmann equation l|19f) _ or by the acceleration equation l|20|) . In order to 
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find the time evolution of the energy density and of the scale factor from the 
system of differential equations one has to prescribe initial conditions for these 
fields and to close the system by choosing an equation of state which relates 
the pressure of the cosmological fluid to the energy density on the brane, i.e, 
p = p(p). 

For the determination of the solution of the system of differential equations 
described above we introduce the dimensionless quantities 

t =~7= t > P=~> P=~, a =— — , (21) 
Vl8 o a a (0) 



and write the conservation equation for the energy density (|14f) , the Friedmann 
(|19ll and the acceleration (|20|l equations in terms of the dimensionless quantities 
as 

p + 3^(p + p) = 0, (22) 
ao 

-J =p+ y + 3 + 4' (23) 

^ = -|(p + 3 P )-^(2p + 3p) + ^-4. (24) 
In equations (|23|l and l|24(l the following abbreviations were introduced 



A = 9 



A p B 



_ 18C 

x ~MoR^- (25) 



The constant A in the Friedmann i|23|) and acceleration (|24l) equations can be 
interpreted as a cosmological constant on the brane. Moreover, it is noteworthy 
to call attention to the fact that equations (|23|) and (|24(l can be written as 

— = P + Prad + + o' ( 26 ) 

ao / 2 3 

— = -^(P + Prad + 3p + 3p rad ) - ^(2p + 3p) + -, (27) 
ao I 16 

thanks to the well-known relationships valid in the four-dimensional case where 
the energy density of the radiation field p ra d scales as 1/oq and the radiation 
pressure is related to the energy density by p ra d = Prad/3. 

From now on we shall analyze the system of differential equations consisting 
of the conservation equation for the energy density \\1'2\\ and of the acceleration 
equation l|24() . In order to close this system of differential equations we have to 
choose an equation of state that relates the pressure to the energy density of the 
cosmological fluid on the brane. In the four-dimensional case one normally uses 
a barotropic equation of state p = wp with < w < 1 to represent a radiation or 
matter dominated Universe. The barotropic equation of state with — 1 < w < 
is also used to represent the scalar fields: the inflaton in the inflationary period 
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of the Universe, and the quintessence in the dark energy dominated epoch of the 
Universe. Other equations of state are also used to model the cosmological fluid, 
namely the Chaplygin equation of state (see e.g. the work [S] and the references 
therein) which can describe a transition from a matter dominated period to a 
cosmological constant dominated epoch and the van der Waals equation of state 
which can simulate the transition from an inflationary period to a matter field 
dominated epoch 0. Here we interested in the description of an accelerated 
period followed by a decelerated epoch, hence we shall use the van der Waals 
equation of state 

P =|^-3p 2 , (28) 
P 

where the parameter w could be identified with the coefficient of proportionality 
in the barotropic formula, since for small values of the energy density p oc wp. 
In classical thermodynamics (see e.g. ^0]) the equation is a reduced van 
der Waals equation of state where the free parameter w is connected with a 
dimensionless temperature. Capozziello and co-workers 0J |S] have used a non- 
reduced form of the van der Waals equation which is characterized by three free 
parameters instead of only one. In their works 0] |S] they also write a reduced 
form of the van der Waals equation with only one free parameter which is not 
similar to the reduced van der Waals equation of state normally found in the 
literature, since the free parameter multiplies also the term proportional to p 2 . 

Apart from the initial conditions for the fields of energy density p(0), scale 
factor a (0) and velocity do(0), one has to specify values for the parameters 
A, x and w m order to find a solution of the system of differential equations 
1|22(1 and Ij24(l which is closed by the van der Waals equation of state (|28|l . 
Here we have specified the initial conditions (by adjusting clocks): p(0) — 1 for 
the energy density, ao(0) = 1 for the scale factor and d(0) = ^/3/2 + A/3 + x 
for the velocity field, which is a consequence of the Friedmann equation (|23() . 
Furthermore, in order to plot the time evolution of the fields in figures 1 and 
2 we have chosen the following values for the parameters A = 0.03, x = 0.5, 
and three different values for w, namely, w = 0.51, w = 0.501, and w — 0.5001. 
Later on we shall discuss how the changes in the parameters w, A and x affect 
the solutions of the system of differential equations. 

In figure 1 it is plotted the acceleration field do (t) as function of time t for 
w = 0.51 (straight line), w = 0.501 (dotted line), and w = 0.5001 (dashed line). 
We infer from this figure that the acceleration field for these three values of w 
describes the distinct accelerated-decelerated-accelerated phases of the Universe. 
In the first period, which refers to a scalar field dominated epoch, the positive 
acceleration grows up to a maximum value followed by a decay towards zero. 
The next phase is related to a matter dominated period where the acceleration is 
always negative and decays to a maximum negative value followed by a growth 
towards zero. The third period is a cosmological constant dominated epoch 
where the acceleration field assumes a positive value. The energy density p 
and pressure p fields are plotted in figure 2 as function of time t for w = 0.51 
(straight line), w — 0.501 (dotted line), and w = 0.5001 (dashed line). We 
conclude from this figure that for the three values of w the energy density field 
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Figure 1: Acceleration d vs time t for w = 0.51 (straight line), for w = 0.501 
(dotted line), and for w = 0.5001 (dashed line). 




Figure 2: Energy density p and pressure p vs time t for w = 0.51 (straight lines), 
for w = 0.501 (dotted lines), and for w = 0.5001 (dashed lines). 



decays with time whereas the pressure field has two distinct behaviors. At the 
beginning the pressure is negative so that the van der Waals fluid behaves like 
a scalar field and it is the responsible for the initial accelerated period. At later 
times the pressure becomes positive so that the van der Waals fluid behaves like 
a matter field and it is partially responsible for the decelerated period, since the 
radiation field which scales as 1 ja% also contributes to the decelerated phase of 
the Universe. When w approaches the value of w = 0.5, we can also infer from 
the figures that the pressure field of the van der Waals fluid at the beginning- 
behaves like an inflaton with an equation of state given by p = — p and it is the 
responsible for an increase of the early acceleration. 

We have chosen three nearby values for w in order to show how the accel- 
eration, energy density and pressure fields behave for small changes of w. We 
proceed now to discuss how significant changes in the values of w have influence 
on the behavior of these fields for fixed values of A and Xi nere ^ — 0.03 and 
X = 0.5. For w = 0.5 the energy density remains constant so that the van der 
Waals fluid behaves like an inflaton with an equation of state given by p = —p. 
the acceleration field grows exponentially and in this case there exists only an 
accelerated phase for the Universe. For values of w < 0.5 the energy density 
grows with time, and we infer that this behavior does not represent a physical 
solution for the evolution of the Universe. For values of 0.5 < w < 0.58 the ini- 
tial accelerated period decreases, since the negative part of pressure of the van 
der Waals fluid decreases, whereas its positive part increases. Hence in the in- 
terval 0.5 < w < 0.58 the three periods accelerated-decelerated-accelerated are 
present in the evolution of the Universe. These three periods for the Universe 
can be found for a larger interval for w by decreasing the amount of radiation, 
i.e., by decreasing the value of the constant x- For values of w > 0.58 (and 
A = 0.03, x = 0.5) the acceleration field evolves from a matter dominated Uni- 
verse where a < to a cosmological constant dominated Universe where a > 0. 
In this last case there exists no inflationary period. 

Let us comment on the behavior of the fields when we change the values 
of the parameters x an( i ^- -By increasing the value of x there exists a more 
pronounced predominance of the radiation field which scales as 1/ap and the 
deceleration period begins at earlier times. By decreasing the value of A, which 
is connected to the cosmological constant on the brane, the present accelerated 
period begins at later times. Moreover, the limit A — > leads to an initial 
accelerated period followed by a decelerated period without a present accelerated 
period, indicating that A plays the role of the dark energy. We could also 
obtain the present accelerated period - which is characterized by a dark energy 
dominated Universe - by following the same methodology of the work [H] and 
instead of introducing a cosmological constant which represents the dark energy, 
we could model the dark energy by an equation of state like the Chaplygin or 
the quintessence equations of state. 

As a final comment we have also investigated the solution of the system of 
differential equations l|22fl and (|24|l closed by a barotropic equation of state for 
the cosmological fluid on the brane, i.e., p — wp with — 1 < w < 1. The results 
we have found for the acceleration field by changing the values of w describe 
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either an accelerated phase, or a decelerated epoch followed by an accelerated 
period of the Universe, i.e., with the barotropic equation of state it was not 
possible to describe the three phases accelerated-decelerated-accelerated found 
here by using the van der Waals equation of state. 
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